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Abstract 

We discuss, using the Hilbert basis method, how to efficiently construct a 
complete basis for L>-flat directions in supersymmetric Abelian and non-Abelian 
gauge theories. We extend the method to discrete {R and non-R) symmetries. 
This facilitates the construction of a basis of all superpotential terms in a theory 
with given symmetries. 



1 Introduction 



Holomorphic gauge invariant monomials play an important role in the understanding 
of supersymmetric theories and phenomenological applications. They represent D-flat 
directions in supersymmetric gauge theories [1] and constitute possible superpotential 
terms. However, in somewhat complex theories the explicit constructions of these mono- 
mials can be quite cumbersome in practice. For instance, already the construction of all 
gauge invariant monomials for the minimal supersymmetric extension of the standard 
model (MSSM) is rather involved j2] (see also [3J for the discussion in stringy extensions 
of the MSSM). 

In [1] it was shown how to construct the so-called Hilbert basis for holomorphic 
gauge invariant monomials which are known to be in one-to-one correspondence 
with the D-fiat directions. The Hilbert basis allows us to write every monomial ^ in 
the form 



Here H is the number of independent basis monomials jfti which can only be determined 
algorithmically. 

The purpose of this note is to extend the notion of the Hilbert basis such as to include 
discrete R and non-i? symmetries as well. In section [2] we start by reviewing the Hilbert 
basis method for continuous gauge symmetries. Section [3] is devoted to the extension to 
discrete symmetries. The general case is discussed in section HJ In section Owe comment 
on potential applications, and finally, section [H] contains our conclusions. 

2 Review of Hilbert bases for continuous gauge sym- 
metries 

Let us briefly review the Hilbert basis method for (continuous) gauge symmetries [I]. 
We start by looking at a theory with a single U(l) factor and then extend the discussion 
to L different U(l) factors. 

2.1 Warm— up example: a single U(l) 

Consider a U(l) gauge theory and fields <jr-f> (1 < / < F) with charges q^'. A monomial 



H 



.ft = \~[ with rji G N . 



(1.1) 




(2.1) 



is gauge invariant if 

q {1) m H h q {F) n F = . 

This condition may be recast as 

rp 

q • n — , 



(2.2b) 



(2.2a) 
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q T = (q^\ . . . , q^ F ') and n T = (rii, . . . , njr). That is, the vector n has to be orthogonal 
to the charge vector q. The requirement that ^# be holomorphic amounts to demanding 
that rii e No- The solutions are the intersection of the hyperplane perpendicular to q 
and the lattice points in Nq 7 . The so-called Hilbert basis 

Jf = {hW,...hW}, (2.3) 

is a complete set of vectors with the property that each solution n of (I2.2p can be 
written as 



H 

n = 

»=i 



Vi h(i) with e N . (2.4) 



Every element of the Hilbert basis is in one-to-one correspondence with a gauge 
invariant monomial 



Jk = (0 (1) ) hil) ---(0 (F) )^ (2.5) 
such that every gauge invariant monomial is given by ( II. ip . 

2.2 Generalization to L U(l) factors 

In the case of L U(l) factors the condition (12. 2p can be rewritten as 

Q-n = 0, neNf (2.6) 
with the charge matrix 

\ 



Q 





■ Qi 


■ 


(F) 



(2.7) 



where q\^ denotes the £ th U(l) charge of the f th field 

An important fact about the above problem is that it is well-known in the mathe- 
matical literature [SHZ] • There are efficient algorithms such as the ones provided by [HIE] , 
allowing us to compute the Hilbert basis for a given matrix Q very efficiently 



2.3 Non— Abelian symmetries 

In the case of a non-Abelian symmetry G, gauge invariance of a monomial composed of 
G representations is equivalent to gauge invariance w.r.t. the r U(l) factors generated 
by the Cartan generators of G. That is, the Hilbert basis method allows us immediately 
to construct G invariant monomials. Some of the resulting monomials are zero and 
others are redundant. Examples for vanishing monomials include the baryons of SU(iV c ) 
theories, e a i a ^-- a N c ^fo) . . . (^* c \ which vanish if, say, i\ = i<i- In order to construct 
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only non-vanishing and inequivalent monomials, there are various methods available. 
One of them is 'consecutive basis building' and the other is to systematically remove 
redundant monomials from the outcome of the Hilbert basis construction. Both methods 
are briefly reviewed in appendix [X] A mathematica package allowing for an automatized 
computation of the independent monomials can be found at [10J. 



3 Discrete symmetries 

In what follows, we generalize the Hilbert basis method such as to applicable to discrete 
symmetries also. This allows us, in particular, to identify a complete basis for allowed 
superpotential terms. 



3.1 Discrete non— R symmetries 

We now discuss how to construct the Hilbert basis for discrete non-f? symmetries. We 
illustrate our method by a simple example, which we then generalize. 



3.1.1 Warm— up example 

We start by a single TLm symmetry under which the cj)^ have charges p^\ The require- 
ment of Ijm invariance of the monomial (12.11) translates into 

mod M . (3.1) 



Hi + 



_|_ p(F) np 



Without loss of generality we can assume that all discrete charges pv> are non-negative. 
Equation (13 .ip is equivalent to 



(-M,pU,...p^) 



( m \ 
\n F J 



with m G N 



(3.2) 



and rif G No, as before. One can now compute the Hilbert basis for the above problem. 
The basis monomials are then given by 



i) 



)' 



■,(0 



'f+i 



(3.3) 



after truncation of the first element hi which represents not a field <f>^ but 



m. 



3.1.2 Multiple discrete non— R symmetries 

The extension to %u\ x • • • x Jjm k is straightforward. We assume that 7jm 1 x • • • x 7jm k 
is already the smallest irreducible symmetry of a possible larger discrete group (cf. [TT]). 
We define the charge matrix 

C = (-M | P) , (3.4) 
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where 



M 



and 



diag(M 1 ,...,M K ) 



Pi 



with p\ denoting the T>M k charge of 
non-negative. 



(3.5) 



(3.6) 



( f) 

We can again assume that all charges p); are 



3.2 Discrete R symmetries 

Let us now turn to discrete R symmetries. We start by discussing a single and then 
generalize the setting to more than one discrete R symmetry. 



3.2.1 Warm— up example 

We start by a single Zjy symmetry under which the (f)^' have charges r(f\ which, again, 
can all be chosen non-negative. The requirement of 1^ invariance of the monomial ( 12 . 1 p 
translates into 



r^ 1 ' Hi + 



2 mod N . 



(3.7) 



Here we have adopted the convention that the superpotential has R charge 2. The above 
equation is equivalent to 



/ £ \ 

m 

m 



with 



(3i 



m G N and n/ G N , as before. One can now compute the solution of the above problem. 
First, one identifies a basis of all vectors orthogonal to (—2, —N,p^\ . . .p^ F ') ■ To deal 
with the restriction £ = 1 we identify two subsets of the Hilbert basis. One, where the 
first entry equals zero, which we call the homogeneous solution space. And one where the 
first entry equals one, which we call the inhomogeneous solution space. If a homogeneous 
solution is labeled by n^ m G M 1 and a inhomogeneous solution by n^ hom G ffi we find 
the general solution to (13. 8p to be 



Ho 



n 



n {i) + V 

'nnhom ~ / j 



(h) 



(3.9) 



h=i 
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with fixed i and rjh G No- That is, there are Hi branches of solutions, where Hi denotes 
the number of inhomogeneous solutions. Each branch consists of Hq solutions with Ho 
denoting the number of homogeneous solutions. We can find all monomials spanning 
the superpotential by truncating the vectors n accordingly. 

To see what that means in practice, consider a setting with a lif symmetry and 



fields ift and (J) with R charges 
example reads 



1 and = 0. The generalized charge matrix for this 



C 



-2,-4,1,0) 



(3.10) 

There is only one vector 
1,0,2,0) T , such that the 
0) T . Similarly, we obtain 



where the last two entries are the R charges of ift and 0. 
orthogonal to C with first component equal to 1, namely ( 
unique inhomogeneous solution is given by n in ho m = (1, 0, 2 
two vectors orthogonal to C with vanishing first component, namely (0,1,4, 0)"* and 
(0, 0, 0, 1) T . The entries of the vectors represent the exponents of the fields in the corre- 
sponding monomials. We hence have found that all allowed superpotential terms are of 
the form 



^2 ^ 2 



(3.11) 



with rji G N . Of course, one could have obtained the result without the Hilbert basis 
method; however, for more complex systems this method is highly advantageous. 

3.2.2 Multiple discrete R symmetries 

Let us extend the discussion to Zj^ x . . . . Define the charge matrix 



(3.12) 





( " 2 




\ 


c = 




-N 


R 




V " 2 




J 



with 



N = diag(iVj 



and 



R 



Here rj^ denotes the j th R charge of the f th field. We compute the kernel of the matrix 
C. The last F components of the kernel vectors with the first entry equal to 1 or 
(before truncation) define the inhomogeneous or homogeneous solutions, respectively. 
The general solution will again be of the form (13. 9p . As before, the identification of the 
(last F) entries of the vectors with the exponents will then give us the desired invariant 
monomials, i.e. allowed superpotential terms. 



N, 



r\ F) \ 



(3.13) 



(3.14) 



5 



4 Putting all together 



4.1 Charge matrix for TJ(1) L symmetry with discrete R and 
non— R symmetries 

We consider now the general U(1) L x Jj Ml x • • • x 1>m k x x • • • ^n, case. The charge 
matrix is 



/ 


-2 






\ 






—N 





R 




-2 














— M 


P 


V 











Q J 



with Q, P, R, N and M defined in equations ([221), M, ( BHD , (B3SD and fl3~5]) . 
respectively. We compute the kernel of C, and, again, we decompose these vectors in 
those with first components equal to 1 or 0. From those we obtain the inhomogeneous or 
homogeneous solutions, respectively, by projecting on the last F components. As before, 
the general solution is of the form (13.91) . and identifying the entries of the vectors with the 
exponents will then give us the desired invariant monomials, i.e. allowed superpotential 
terms. 



4.2 A stringy example 

We base our example on the Z 2 x Z 2 orbifold discussed in [4] . We consider the Zf vacuum 
discussed there and construct the superpotential for the standard model singlets with R 
charges or 2. At the orbifold point, the symmetry seen by these fields is of the type 

Gsymm = U(1) L X (Z 2 ) 6 X (Zf ) 3 (4.2) 

with L = 8. Here we have already eliminated the non-Abelian symmetries by forming 
SU(iV) gauge invariant monomials (cf. our discussion in 12. 3p . The discrete symmetries 
follow from the space-group selection rules and i^-momentum conservation [12J (for the 
rules in this specific geometry see [13]), and are partially redundant, e.g. there are in 
fact only two independent lif factors. These symmetries and selection rules constrain 
the holomorphic correlators of the theory [T2"tlT4"] , which can also partially come from the 
Kahler potential in an appropriate description (cf. the discussion in [T5J section 4]). In 
what follows we will not distinguish between such correlators and allowed superpotential 
terms. 

In the vacuum discussed in |4], there is a residual lif symmetry, which forbids the 
superpotential at the perturbative level. If we switch on an additional field with R 
charge 2 we will break the and obtain a non-zero superpotential in the vacuum. We 
will switch on the additional field <p\. That is, the fields 

4> (l) = {01, 01, 02>3, 04, 05, 06>7, 08, 09>1O, 011, 012, 013, 014, ^ ^ 

xi,x 2 , x 3 , x 4 , x 5 ,xi, x 3 , x 4 , x 5 , y 3 , y 4 , J/5, He} 
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will now acquire vacuum expectation values (VEVs). The Hilbert basis contains 15408 
elements for this choice. The superpotential starts at lowest order with 4 Hilbert basis 
elements 

W = (£4X4 + X 5 X 5 + 9 0i3 + 0io01 4 ) 01 H • (4.4) 

We also considered the appearance of the proton decay operator QQQi. The Hilbert 
basis involving QQQi consists of 4284 elements where we focus only on first generation 
quarks and leptons. The lowest order QQQi operator occurs at order 11 in the field 
VEVs. An example is given by 

W D Q1Q2 Q2 t\ 01 Xl X 2 X 3 X 4 X 3 02 04 09 0^2 ■ (4-5) 

These examples show that the Hilbert basis method is powerful enough to handle very 
complex examples with many fields. 



5 Applications and speculations 



Defining the subsets J^ mhom = {^Si! om , • • • , ^ in hom} and = {-^w ■ ■ ■ > -^S) 

^inhom n hom' 

Hi 



which are constructed from n^L and n^' ' , the full superpotential to all orders is 



i ni 



We can further speculate that the structure of the superpotential is 



Hi m 

^inhom 

•^structure • (^'^) 

h=i 

However, the relations between the couplings implied by W stTUCture will in general be 
incorrect. Yet Structure can be used quickly to answer certain questions such as at which 
order some combination of fields appears first. Such questions arise in (generalized) 
Froggatt-Nielsen model building [16J with many symmetries and where the order in 
which a given term appears is a measure for the suppression of the corresponding effective 
coupling. Another, very similar application concerns the identification of approximate 
continuous R and non-i? symmetries, where one is interested in the lowest order terms 
that explicitly break such symmetries [T7TfT8] . Given (I5.2p one immediately reads off the 
lowest order at which perturbative superpotential terms arise. 

Whether or not the relations between the coefficients implied by (15. 2p can be ob- 
tained in some interesting setting needs still to be explored. However, it is tempting to 
speculate that in certain highly symmetric string compactifications, such as orbifolds, 
the superpotential (and holomorphic terms in the Kahler potential) will be some known 
function of the building blocks. It will be interesting to study this question in more 
detail. 



7 



6 Conclusions 



We have described a simple method that allows us, given the symmetries of the theory, 
to construct the building blocks of the .D-flat directions and the superpotential. This 
basis is given by the basis of non-negative integer solutions n of the simple matrix 
equation C ■ n — 0, which has been extensively studied in the mathematical literature. 
Publicly available codes allow us then to compute the basis very efficiently We have 
discussed a specific example, based on a Z 2 x Z 2 orbifold model, in which we computed 
the superpotential basis for certain singlet fields. The successful construction of this basis 
demonstrates that our methods can be used to efficiently compute the superpotential to 
all orders even in rather complex systems. 
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A Gauge invariant monomials for SU(AT) 

We review the construction of gauge invariant monomials for SU(iV) with matter fields 
(pi and <pj. 

A.l Consecutive basis building 

Let us consider a general example of L SU(iVj) gauge groups, SU(A r 1 ) x • • • x SXJ(N L ). 
In order to construct the Hilbert basis M' we proceed as follows: in a first step, we 
construct a basis M{ of SU(iVi) singlets, consisting of elementary SU(iVi) singlets and 
SU(iVi) invariant monomials. As is well known, the latter will be given by the 'mesons' 
and 'baryons', which in SU(3) would look like, 

These monomials will transform as singlets under SU(A^!) and together with the singlet 
fields they build the basis M{. 

We now use J^i to build the basis which will be the basis of SU(A r 1 ) x SU(iV2) 
singlets. Obviously, will contain terms of which are already SU(iV2) singlets and 
also monomials which are constructed in a similar way as in (lA.lj) . The only difference 
is that terms in M{ can have more than one SU(iV2) index, i.e. they can furnish higher 
represent at ions . 



S 



From here on the course of action is always the same: we use the previous basis 
3ti[,...,k to build J0|,„.,fc,k + i, the basis of SU(iVi) x • ■ • x SU(iV fc ) x SU(iV fc +i) singlets until 
we eventually find the basis Ji? = <ffi\ l of monomials invariant under the full gauge 
group. One can find an explicit example of this method for SU(3) x SU(2) x U(l) in [2J. 

A. 2 Cartan subalgebras 

Alternatively, one can construct the monomials for non-Abelian symmetry groups by 
using the Hilbert basis method for the U(1) L symmetry defined by the Cartan subalge- 
bras. For this strategy, we split all fields in their tensor components and assign them a 
charge according to the Cartan charges. Let us consider some fields from [2] to explain 
this procedure. Assume, we have an SU(3) x SU(2) gauge group. SU(3) has rank 2 
and its Cartan subalgebra can be taken to be the one generated by the two diagonal 
Gell-Mann matrices, A3 and Ag, whereas SU(2) has rank 1, hence we use the diagonal 
Pauli matrix 03. 

Let us assume we have three fields, Q(3, 2), u(3, 1) and £(1,2). Q has six tensor 
components Q m , where a = 1, 2, 3 is the SU(3) index and i = 1, 2 the SU(2) index. Now 
we can assign charges to each component under the respective gauge groups, using the 
eigenvalues of the generators of the Cartan subalgebras, which is particularly easy when 
using the diagonal matrices A3, A 8 and a 3 . Therefore, Q 2 = Q a=2,l=1 will be assigned the 
charges 

= - 1 , qf = 1 , q^ = 1 . (A.2) 

Another example is u a , which is an SU(2) singlet and an SU(3) anti-triplet, thus 
carrying zero charge under the SU(2) gauge group and opposite SU(3) charges. Taking 
the component TL2 = u a=2 , we get 

?? 2 = 1, q? = ~l, qf = 0. (A.3) 

In this way we can split each field into its components (six for Q, three for u and two 
for £) and build a 3 x 11 charge matrix Q, which in our example will look like 

/ 1 -1 1 -1 -1 1 \ 
Q = 1 1 -2 1 1 -2 -1 -1 2 0. (A.4) 
\ 1 1 1 -1 -1 -1 1 -1 / 

Using algorithms like [SIE], we can now find the Hilbert basis of all solutions for the 
charge matrix Q, the same way as before in equation ( 12. 7p . Solutions would for example 
be Q n Q 21 Q 32 £ 2 , u 1 !! 2 ^, Q 11 !! 1 ^ or many more. 

In order to associate these solutions to proper gauge invariant monomials, we define 
a prescription of how to translate such expressions to objects in which the (generalized) 
color indices are contracted appropriately. An increasing series will be contracted with 
the total antisymmetric Levi-Civita tensor, indices which have the same value but are 
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dotted and undotted will be contracted with the Kronecker delta. Using this procedure, 
our examples would look like 



Q u Q 21 Q 32 ( 2 ^ Q aa Q pb Q? c i d e aPl e ac e bd , (A.5a) 
where we had to contract antisymmetrically several times, 

rtlfui xPvPute^, (A.5b) 
which is very similar to the one above and 

Q ll uH 2 ^ Q aa v*£ b 5 adl e ab , (A.5c) 

where we were able to see when to use the Kronecker delta. 

A caveat of this procedure is that we will end up with many monomials occurring 
more than once, e.g. flA.5bj) will appear six times. Therefore, we have to remove the 
redundant ones. Furthermore, one has the possibility to end up with monomials which 
will vanish due to antisymmetric contraction. Take flA.5bj) again: if u would be a field 



with only one generation, the monomial would clearly vanish. This means, one has to 
check the Hilbert basis for zero-valued monomials and remove them, which will leave 
the basis nonetheless intact, since the contribution of these monomials would be zero in 
any case. 

We see that all these caveats are manageable. Furthermore, this procedure has a 
big advantage compared to the method described in the previous subsection IA.1I Using 
the Cartan subalgebras allows us to quickly and fully automated build monomials for 
several SU(iV) and U(l) gauge groups, which would get very tedious (in certain cases 
impossible), especially for N > 3, more than three gauge groups or too many fields. We 
have created a mathematica package ourselves, based on [SJIS], which automatizes this 
procedure [TU] . 
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